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ABSTRACT 

We present the X-ray point source population in the nearby Virgo elliptical galaxy NGC 4636 from three 
Chandra X-ray observations. These observations, totaling ~ 193 ks after time filtering, were taken with the 
Advanced CCD Imaging Camera (ACIS) over a three year period. Using a wavelet decomposition detection 
algorithm, we detect 318 individual point sources. For our analysis, we use a subset of the 277 detections with 
> net 10 cts (a limiting luminosity of approximately 1.2 x 10 37 erg s" 1 in the 0.5 - 2 keV band, outside the 
central 1 .5' bright galaxy core). We present a radial distribution of the point sources. Between 1 .5' and 6' from 
the center, 25% of our sources are likely to be background sources (AGN) and 75% are LMXBs within the 
galaxy, while at radial distances greater than 6', background sources (AGN) will dominate the point sources. 
We explore short and long-term variability (over time scales of 1 day to 3 years) for X-ray point sources in 
this elliptical galaxy. Fifty-four sources (24%) in the common ACIS fields of view show significant variability 
between observations. Of these, 37 are detected with at least 10 net counts in only one observation and thus 
may be "transient". In addition, ~10% of sources in each observation show significant short term variability; 
we present an example light curve for a variable bright source. The cumulative luminosity function for the 
point sources in NGC 4636 can be represented as a power-law of slope a = 1.14 ±0.03. We do not detect, 
but estimate an upper limit of ~ 4.5 x 10 37 erg s" 1 to the current X-ray luminosity of the historical supernova 
SN1939A. We find 77 matches between X-ray point sources and globular cluster (GC) candidates found in 
deep optical images of NGC 4636. In the annulus from 1 .5' to 6' of the galaxy center, 48 of the 129 X-ray point 
sources (37%) with > 10 net counts are matched with GC candidates. Since we expect 25% of these sources to 
be background AGN, the percentage matched with GCs could be as high as 50%. Of these matched sources, we 
find that ~ 70% are associated with the redder GC candidates, those that are thought to have near-solar metal 
abundance. The fraction of GC candidates with an X-ray point source match decreases with decreasing GC 
luminosity. We do not find a correlation between the X-ray luminosities of the matched point sources and the 
luminosity or color of the host GC candidates. The luminosity functions of the X-ray point sources matched 
with GCs and those that are unmatched have similar slopes over 1.8 x 10 37 erg s -1 < L x < 1 x 10 38 erg s -1 . 
Finally, we present a color-color diagram based on ratios of X-ray flux rather than source counts, which yields 
a much tighter source distribution, and shows a large population of sources which are likely LMXBs and a 
small population of black hole candidates. 

Subject headings: galaxies: individual (NGC 4636) — X-rays: galaxies — X-rays: binaries 



1. INTRODUCTION 

The sub-arcsecond resolution of the Chandra Observatory 
has revealed the nature of X-ray emitting point sources in 
nearby galaxies. It is clear that almost all the very lumi- 
nous {Lx > 10 36 erg s" 1 ) point sources in galaxies belong to 
two distinct populations of compact binaries, their evolution- 
ary timescales depending upon that of their donor stars: the 
low-mass X-ray binaries (LMXB), which are long-lived and 
evolve on timescales of 10 9 - 10 10 yr, and high-mass X-ray bi- 
naries (HMXB), which evolve on timescales of 10 6 - 10 7 yr. 
The latter population is thus an indicator of recent star forma- 
tion and is not expected to be found in early-type galaxies, un- 
less a recent merger has occurred. The LMXB population, on 
the other hand, has lifetimes comparable with that of the host 
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galaxy, and their number and combined luminosity is found 
to correlate well with the stellar mass of galaxies (Gilfanov 
2004). 

In this paper, we present a Chandra view of the X-ray point 
source population of NGC 4636, a bright E/SO galaxy on the 
southern periphery of the Virgo cluster. It has a radial veloc- 
ity similar to that of Virgo, but is 10.8° from the center of the 
cluster, which corresponds to 2.8 Mpc at a distance of 15 Mpc 
(Tonry et al. 2001). Furthermore, the galaxy lies at the cen- 
ter of a poor group (Osmond & Ponman 2004; Miles et al. 
2004, 2005; Baldi et al. 2008), possibly falling into the clus- 
ter. Its unusual properties have attracted detailed multiwave- 
length research for several decades. It has been suggested 
that the galaxy has an unusually large dark halo (Loewenstein 
& Mushotzky 2003, Schuberth et al. 2006, Chakrabarty & 
Raychaudhury 2008). NGC 4636 was one of the first early- 
type galaxies in which neutral hydrogen was detected (Knapp 
et al. 1978), and further radio observations (Birkinshaw & 
Davies 1985) revealed a weak central radio source and small- 
scale jets. NGC 4636 has a large population of globular clus- 
ters similar to ellipticals of comparable luminosity (Dirsch et 
al. 2005). Its far-infrared luminosity greatly exceeds that ex- 
pected from its stellar content (Temi et al. 2003), and the lu- 
minosity in the vicinity of 100/i is consistent with dust emis- 
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TABLE 1 

Summary of Chandra/ ACIS Observations of NGC 4636 



optical DSS image in FigureQ] 



Date 


Detector 


Sequence # 


Exposure Time (s) 


2000 Jan 26 


ACIS-S 


600083 


44450 


2003 Feb 14 


ACIS-I 


600300 


74709 


2003 Feb 15 


ACIS-I 


600331 


74190 



sion from a recently accreted disk galaxy. Its flattened (E4) 
morphology at the outer faint isophotes (Sandage 1961) in- 
dicates the presence of large-scale angular momentum, often 
associated with recent mergers. 

NGC 4636 is one the brightest nearby early-type galaxies 
in X-rays, and so it has been well-studied with generations of 
X-ray observatories. It was first detected as an extended X-ray 
source by Einstein (Forman et al. 1985). ROSAT and ASCA 
observations found abundance and temperature gradients in 
its extended X-ray halo (Trine hieri et al. 1994; Matsushita 
et al. 1997, Buote 2000), while Chandra and XMM-Newton 
observations (Jones et al. 2002; O' Sullivan et al. 2005) show 
symmetric arm features and cavities in the extended X-ray 
halo, interpreted to be evidence of past AGN activity. 

The resolution of the Chandra observatory has enabled the 
detailed study of the environment of detected point sources 
(e.g. Fabbiano 2006 and references therein). Many of these 
are associated with globular clusters. Since it has been sug- 
gested that LMXBs may be primarily formed in the cores of 
globular clusters (White et al. 2002), the properties of point 
sources, along with those of their host GCs, may yield im- 
portant clues to the role played by GCs in the formation of 
LMXBs (e.g. Maccarone et al. 2004. Kundu et al. 2007). On 
the other hand, there is some evidence that a significant frac- 
tion of LMXBs may be formed in the field, and are thus not 
expected to be associated with GCs (Juett 2005, Irwin 2005). 
The collective study of X-ray detected point sources and glob- 
ular cluster candidates found in optical studies can thus yield 
important information concerning the origin of X-ray emitting 
binaries in galaxies. 

We introduce the observations and d iscu ss data reduction 
and point source detection in $2] In jj3.ll we examine the 
distribution of sources in NGC 4636. In 83.21 we use soft, 
medium, and hard band source fluxes to make an X-ray color- 
color diagram whic h is independent of any instrumental or 
detector effects. In 83.31 we examine the X-ray point source 
population to l ook f or variable sources on both long and short 
timescales. In 83.41 we pr esent the luminosity function of X- 
ray point sources, and in 83.51 we optically identify sources 
associat ed w ith globular clusters and examine their proper- 
ties. In 83.61 we examine a few interesting individual sources. 
Finally, our analysis and results are summarized in §0] 

2. OBSERVATIONS, DATA PROCESSING, AND POINT 
SOURCE DETECTION 

The observations presented in this paper were made with 
the ACIS-I and ACIS-S detectors on Chandra on three occa- 
sions spanning three years. The observations are summarized 
in Table Q] 4 The FOV of each observation is overlaid on an 

4 We also analyzed a short (5 ks) ACIS-I observation (sequence #600084) 
taken in December 1999, but it yielded no sources that were not detected in 
the other observations, and the source counts in the few detections present 
were too low to allow meaningful analysis. Thus, we do not use this observa- 
tion in the work discussed in the remainder of the paper. 




FIG. 1. — A DSS optical image of NCG 4636 with the FOV of our Chandra 
observations overlaid. The light blue outline is the January 2000 ACIS-S ob- 
servation, and the green outlines are the February 2003 ACIS-I observations. 



We obtained level 1 event lists for the observations from 
the Chandra data archive. We filtered these event lists to in- 
clude only the standard good grade set (ASCA grades 0,2,3,4 
and 6) and to exclude bad pixels. The event lists were also 
filtered on the nominal good time intervals for each observa- 
tion. Using the CIAO Contributed S-Lang Script lc_clean.sl, 
we did additional time filtering to remove times affected by 
background flaring. We set the script so as to first calculate 
an initial mean (3 sigma-clipped), then consider times with a 
count rate greater than 1.2 times this mean to be flaring. The 
remaining total good exposure times are listed in Table Q] 

We created a 0.5 - 2 keV band image for each observa- 
tion and used the wavelet decomposition algorithm wvdecomp 
(Vikhlinin et al. 1995) on the images to detect sources above a 
threshold of 4.5er. With this detection threshold, and given the 
image sizes, we expect approximately 2 false detections per 
image (Vikhlinin et al. 1995). We considered a range of scales 
from 1 to 6 image pixels (equivalent to roughly 1 to 6 arcsec- 
onds, since 1 ACIS pixel ks 0.5" and we binned by a factor 
of 2). In addition, we co-added the two ACIS-I observations 
done sequentially in February 2003, and ran the source detec- 
tion algorithm on this image. After detection, source locations 
were refined using a centroiding algorithm, and the 90% en- 
circled energy radius, based on the preflight calibration model 
of the Chandra PSF, was determined for each source location. 
We manually examined the detected sources, looking at them 
in the ACIS images and in the individual scale wavelet de- 
composition images, and rejected some detections along the 
galaxy arms which appeared to be knots of hot gas rather than 
stellar point sources. We also rejected spurious detections 
along chip gaps. 

For each observation, we constructed a 0.5 - 2 keV band 
background image by summing the large-scale (5 and 6) 
wavelet decompositions with the residual image (i.e. the im- 
age minus the total output of all scales used). For each source, 
we used a background region centered on the source coordi- 
nates with a radius of 1 .5 times the 90% encircled count frac- 
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tion radius for the source. We also made an exposure map for 
each observation. For each source detected in each observa- 
tion, we tabulated raw and net (i.e. background-subtracted) 
counts, and raw and net exposure-corrected counts. We com- 
pared source positions between the observations, and matched 
sources whose coordinates agree within a dynamic matching 
radius dependent on the PSFs at the locations of the sources. 
This dynamic matching radius was set to be equal to the aver- 
age of the 90% encircled energy radii of the two sources, or, 
if this quantity was smaller than 1 .5 arcseconds, the matching 
radius was set to 1.5 arcseconds. 

As a final step in the initial analysis, we calculated a counts- 
to-flux conversion factor for each source by computing re- 
sponse matrices at each source location on the appropriate 
ACIS chip, assuming a power-law model with T = 1 .5 (Irwin 
et al. 2003) and n H = 1.81 X 10 20 cm" 2 (Dickey & Lockman 
1999). After comparing fluxes computed with conversion fac- 
tors based on different values of T and n#, we conclude that 
the choice of parameter values, within a range of V = 0.5-2.5 
and n H = 9.0 x 10 19 -3.6 x 10 20 cm" 2 , changes the resulting 
fluxes by less than 25%. To convert source fluxes to luminosi- 
ties, we assume a distance of d = 15 Mpc (Tonry et al. 2001). 
To facilitate comparison with other studies, we note that the 
counts-to-flux conversion factors for the 0.3-8 keV band are 
typically 2.2 times greater than those for the 0.5-2 keV band 
for the ACIS-S, and typically 2.0 times greater for the ACIS-I. 

From the 0.5 - 2 keV background images, we estimate the 
luminosity a source would need to be detected at a given lo- 
cation. By looking at the maximum values of these "S„„„" 
maps, we estimate that outside of the bright 1.5' central core, 
our source detection is complete above ~ 1.8 X 10 37 erg s" 1 . 
Of the 277 sources that we detect with > 10 net counts, 244 
are above this luminosity. 

A summary of X-ray point sources detected in the NGC 
4636 observations with at least 10 net counts is given in Table 

El 



shape, with the deeper ACIS-I observations having a larger 
number of sources with greater than 10 net counts. 




FIG. 2. — Smoothed, exposure-corrected images in the 0.5 - 2 keV band. 
North is up and East is to the left. The left image is from the co-added ACIS- 
I observations taken in February 2003, totaling approximately 150 ks. The 
right image is from the 45 ks ACIS-S observation taken in January 2000. 
White circles show sources detected in the ACIS-I image, red circles show 
sources detected in the ACIS-S image, and blue circles show sources detected 
in both images. The logarithmic scale (identical in both images) ranges from 
0-1 X 10~ 3 cts sec -1 pixel -1 . 



3. RESULTS 

3.1. The radial distribution of the point sources 

The distribution of the X-ray point sources as a function of 
distance from the center of NGC 4636 is shown in Figure [3] 
The radial profiles of sources in the ACIS-S and ACIS-I obser- 
vations are shown separately, and are seen to have the same 



ACIS-S (600083) 
ACIS-I (600300+600331) 
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Distance from Center (arcmin) 

FIG. 3. — Radial source density profile for sources with > 10 net counts. 
The ACIS-S observation (sequence # 600083) is shown by a dashed line and 
the sequential co-added ACIS-I observations (600300 and 600331) are shown 
with a dotted line. From results of the Chandra Deep Fields (Bauer et al. 
2004), we estimate an AGN density of ~ 0.3 arcmin -2 at our limiting flux of 



'4x10 16 erg cm 2 s 1 



in the 0.5 - 2 keV band. This density is shown by 
the solid line. Note that our field becomes AGN-dominated outside of ~ 6' . 



From point source statistics in the Chandra Deep Fields 
(Bauer et al. 2004), we estimate a background AGN den- 
sity of ~ 0.3 arcmin" 2 at our limiting flux of ~ 4 x 10" 16 erg 
cm" 2 s" 1 in the 0.5 - 2 keV band (shown as a horizontal line 
in Figure[3]). Based on this estimate, our field becomes AGN- 
dominated beyond ~ 6' of the galaxy center. Between 1.5' 
and 6', we estimate that 25% of our sources are background 
sources (AGN) and 75% are LMXBs within the galaxy. 

3.2. Spectral Analysis 

Traditionally, X-ray colors for sources are calculated as ra- 
tios of counts in different energy bands (e.g. Swartz et al. 
2004; Prestwich et al. 2003). However, since we want to di- 
rectly compare sources observed by X-ray detectors with dif- 
ferent responses (i.e. the ACIS-S and ACIS-I) and on different 
parts of the detectors, we convert source counts to fluxes be- 
fore computing X-ray colors. We define three bands: soft (S 
= 0.5-1 keV), medium (M = 1 - 2 keV), and hard (H = 2 
- 8 keV), and calculate two colors as Ci = (M-S)/(S+M+H) 
and C 2 = (H-M)/(S+M+H), as defined in Swartz et al. (2004). 
The resulting color-color diagram is shown in Figure [4] The 
plotting symbol size is proportional to the 0.5 - 2 keV band lu- 
minosity. We identify two distinct populations: a large cluster 
of harder, less luminous sources with power-law indices be- 
tween 1 and 2, and a smaller group of softer, more luminous 
sources (Lx ~ 10 38 erg s" 1 ) with steeper power-law indices. 
Points in black lie less than 6' from the center of the galaxy 
and are most likely members, while points in red lie further 
than 6' from the center. The source represented by a purple 
X is located at the galactic nucleus, and the two sources with 
steep spectra represented by the purple asterisk and plus sign 
lie within 5 " of th e nuc leus. T hese may be signatures of black 
holes (see 33. 6. H and 33.6.2b . The large red circle near these 
points is a bursting X-ray source, which we discuss further in 
33.6.41 Due to its distance from the galactic center (> 1 1') it 
is most likely not a member of NGC 4636. 
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FIG. 4. — Color-color diagram of all point sources in the NGC 4636 ob- 
servations having at least 5 net counts in each of the three bands. We follow 
the example of Swartz et al. (2004) and use bands S = 0.5 - 1 keV, M = 1 - 
2 keV, and H = 2 - 8 keV (T = S+M+H, i.e. 0.5 - 8 keV). Note, however, 
that these color ratios are calculated from source flux, as opposed to counts. 
This allows us to compare ACIS-I and ACIS-S directly. The different plotting 
symbols indicate whether the source was detected in the ACIS-I observation 
(open square), the ACIS-S observation (open diamond) or both observations 
(open circle). The three purple points denote sources detecte d within 5" of 
the galactic center: the X is at the galaxy nuc leus (s ee 4)3.6. H and the + and 
* are two nearby blackhole candidates (see jB.6.21 . Plotting symbol size is 
proportional to source luminosity in the 0.5 - 2 keV band. Sources in red lie 
further than 6' from the center of the g alaxy. The red point near the asterisk is 
a bursting X-ray source, discussed in 33.6.41 Blue curves denote colors of ab- 
sorbed power-law models of spectral indices V = 1, 2, 3, and 4 over the range 
of absorbing column densities rifj = 10 20 to 10 24 cm 4 . Green curves denote 
constant absorption column densities, as given in the legend. A typical error 
bar is shown in the upper right corner. 



To make this flux based color-color diagram, we performed 
an iterative flux calculation as follows. Beginning with source 
counts in each band for both the ACIS-I and ACIS-S de- 
tectors, we convert to source and model flux, assuming a 
power-law model with slope T = 1.5 and column density 
n# = 1.81 x 10 20 cm" 2 and computing response matrices at 
each location on the appropriate ACIS chip (as described in 
§J2}- We convert the model curves from counts to flux as- 
suming the same power-law model and using response ma- 
trices calculated at the aimpoint location of each observation. 
Since in earlier works the color-color diagram has been based 
on counts, we show in Figure [5] the counts-based diagram 
(with model curves for the ACIS-I). We show in Figure |6]the 
color-color diagram from this first flux iteration. Next, we 
use model grids and source colors from this first iteration to 
determine the values of T and n# closest to each source in 
color-color space. We then use these individual values of T 
and Tin to calculate a more accurate counts-to-flux conversion 
factor for the given source (again computing response matri- 
ces at each location on the appropriate ACIS chip). Note that 
the model grid shown in Figure|6]is a subset of the much finer 
grid that we use to pick the individual T = 1 .5 and n# for each 
source. For example, our grid has T ranging from 1 to 4 in 
steps of 0.1, but for clarity on the plot, we only show integer 
values of T. 

After performing this iterative flux conversion for the 
sources, we again calculate the colors, now comparing them 
to grids directly calculated from model flux values. This is 
shown in Figure |4] We note that this figure has much less 
scatter than the color-color plot based on counts for the same 
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FIG. 5. — As in Figure|4] only here the colors and model curves are cal- 
culated from counts, whereas in Figure f4] they are calculated from fluxes. 
The model curves in this plot are based on the ACIS-I response. Note that 
here, unlike in Figuref4] it is not possible to easily group the sources into two 
populations. 
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FIG. 6. — This plot shows the middle step in our iterative flux calculation 
process. Here, the source and model counts in each band have been con- 
verted from counts to flux assuming a power-law model with T = 1.5 and 
n# = 1.81 X 10 20 cirT 2 . This plot shows model curves for the ACIS-I, and 
we have plotted only sources detected in the co-added ACIS-I observations. 
Based on the location of each source in color-color space, we select the near- 
est model curves, and use that V and n# to recalculate the counts-to-flux con- 
version factor for that source. (An identical process is followed for sources 
detected in the ACIS-S observation, using model curves for the ACIS-S.) The 
results of this iterative flux calculation are shown in Figuref4] 

set of sources (Figure [5}, making it easier to identify popula- 
tions and trends. 

3.3. Temporal Variability 

We present a comparison of source fluxes between the Jan- 
uary 2000 ACIS-S observation and co-added February 2003 
ACIS-I observations in Figure|7] To search for long-term vari- 
able sources, we determine a significance threshold as 



S = 



l/s-Z/l 



>3 



(1) 



af + at 



where / is the 0.5 - 2 keV band flux and a is the flux un- 
certainty (dominated by the Poisson error on the number of 
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FIG. 7. — A comparison of source fluxes in the 0.5-2 keV band from the 
45 ks ACIS-S observation and co-added February 2003 ACIS-I observations. 
Only sources with > 10 net counts in both observations are plotted here. 
Variable sources are shown in red. 

counts - the error on the counts-to-flux conversion factor is 
minimal and we do not include it here). For the 228 sources 
in the common field of view detected with > 10 net counts in 
at least one of the observations, we find 54 (~24%) are long- 
term variable sources. These are marked with a "V" in the 
last column of Table|2] Of these 54 sources, 17 have > 10 net 
counts in both observations. These are shown in red in Figure 
[7] The remaining 37 variable sources are "transient" - that is, 
they have > 10 net counts in only one observation, and are 
not "reliably" detected in the other observation. Of these 37 
transient sources, 3 1 have > 20 net counts in one observation 
and < 10 net counts in the other observation. 

To measure variability on shorter timescales, we compare 
source fluxes between the two February 2003 ACIS-I obser- 
vations (see Figure[8]). These observations were done succes- 
sively, each lasting for ^75 ks (~21 hrs). Of 188 sources de- 
tected in the common field of view with at least 10 net counts 
in one observation or the other, 9 (5%) vary significantly (S > 
3) between the two observations. These sources are marked 
"MV" in the last column of Table [2] Figure [8] shows a subset 
of these sources: those detected with at least 10 net counts in 
both observations. 

For yet another look at short term variability, we use 
the IDL PINTofALE program timevarvk (Kashyap & Drake 
2000), to search for sources which vary within the course 
of an observation. The timevarvk routine uses a one sample 
Kolmogorov-Smirnov test to compare the source lightcurve 
to a flat model, then calibrates the observed deviation with 
Monte Carlo simulations. Good time intervals are accounted 
for and gaps are removed. In addition, the deviations between 
the source lightcurve and model are averaged over a given 
number of photons to suppress random Poisson deviations. 
We ran 1000 Monte Carlo simulations per source, and aver- 
aged over 2 photons. Of 132 sources detected in the ACIS-S 
observation with at least 10 net counts, 17 (13%) are variable 
at a 3-sigma level (p < 0.0027). Of 233 sources detected in 
the summed ACIS-I observation with at least 10 net counts, 
19 (8%) are variable at a 3-sigma level. These sources are 
marked with "SV-I" and/or "SV-S" in the last column of Ta- 
ble|2] 

An example of a variable lightcurve is shown in Figure 
[9] This source shows short term variability during both the 
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FIG. 8. — A comparison of source fluxes in the 0.5 - 2 keV band from the 
2003 February 14 75 ks ACIS-I observation (sequence number 600300) and 
the 2003 February 15 75 ks ACIS-I observation (600331). Only sources with 
> 10 net counts in both observations are plotted here. Variable sources are 
shown in red. 



ACIS-I and ACIS-S observations, and also shows variability 
between the two observations. Note the flare in the ACIS-S 
observation (top panel). We estimate its duration to be ap- 
proxim ately 6 00 s. We discuss this bursting X-ray source fur- 
ther in 

Comparing our short and long term variability results to X- 
ray point source surveys of other galaxies is challenging due 
to the different methods and criteria used to determine and 
quantify variability in each case. To minimize the effects of 
these differences, we attempt to compare results for a subset 
of the brightest sources in each galaxy. For example, Kraft et 
al. (2000) find 35 out of 246 (14%) of X-ray point sources in 
Cen A to be variable with a > 3er significance over a 5 month 
period separating the two observations. Computing the chi- 
square statistic for individual light curves (binsize ~ 3600 s), 
they find only 2 sources with short term (within an observa- 
tion - each about 36 ks) variability. Of the 17 sources with 
luminosity greater than 10 38 erg s" 1 , 8 (47%) show significant 
long term variability. Jordan et al. (2004) search for long term 
variability among a subset of sources detected in two M87 ob- 
servations spanning two years. Of the 23 sources with lumi- 
nosity greater than 10 38 erg s , 7 (30%) vary significantly. 
In NGC 4636, we find that of the 52 sources with luminosity 
greater than 10 38 erg s , 21 (42%) show some type of vari- 
ability: 13 (24%) long term and 12 (23%) short term. (Four 
of these very luminous sources show both long and short term 
variability.) 

Loewenstein et al. (2005) find 5 transient sources out of 39 
(13%) and an additional 6 (15%) highly variable sources in 
two observations of NGC 1399 spanning three years. Three 
of these 6 sources have luminosities greater than 10 39 erg s" 1 , 
although since luminosities are not given for the whole source 
list, we have no basis for comparison with our results. 

Of 126 sources with luminosity greater than 1.4 x 10 37 
erg s" 1 in NGC 4697, Sivakoff et al. (2008) find that 5 sources 
(4%) display short term variabilty and 16% display long term 
variability. Restricting our NGC 4636 population to the 157 
sources lying less than 6' from the galactic center with lumi- 
nosity greater than 1.4 x 10 37 erg s" 1 , we find 16 (10%) that 
display short term variability and 29 (18%) that show long 
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FlG. 9. — Variable source detected at a(72000) = 190.8825°, 5(72000) = 2.6190°. This source shows both short term (in both observations) and long term 
(between observations) variability. Counts are in the 0.5 - 2 keV band. Background counts are plotted with a dotted line. The upper panel is from the ACIS-S 
observa tion (sequence # 600083) and the lower panel is from the co-added ACIS-I observations (sequence #s 600300 and 600331). This source is discussed in 
jl3A4l 



term variability. 



3.4. Luminosity function of point sources 

Figure [TOl shows the luminosity function in the 0.5 - 2 keV 
band for X-ray point sources lying between 1.5' and 6' from 
the center of NGC 4636. We plot here only the sources that 
have 10 net counts, after subtraction of the local background, 
in each observation. Our requirement of having at least 10 
net counts per source, leads to a different minimum luminos- 
ity for the two sets of observations (1.2 x 10 37 erg s" 1 and 
1 .8 x 10 37 erg s _1 for ACIS-I and ACIS-S respectively). At the 
bright end, consistent with Irwin, Bregman & Athey (2004) 
and Raychaudhury et al. (2008), we do not find any ULXs 
(L x > lxlO^ergs" 1 ). 

We model the luminosity function of these sources as a sin- 
gle power-law, since we see no evidence for a break in the 
luminosity function. Jordan et al. (2004) demonstrate that 
there is no compelling evidence for a break in the luminosity 
functions of the Virgo cluster galaxies M87, M49, and NGC 
4697. We adopt a robust method of measuring the slope of the 
luminosity function, which is usually expressed in the cumu- 
lative form 



S is 



logN(> S) = -alogS+K. 



(2) 



Each of our n point sources has a measured luminosity F;, 
given the adopted distance to the galaxy, and an estimated 
error a,. The probability of a point source to have a luminosity 



On comparison with d2j, /3 = l + a 
likelihood function 



P(S)dS = AS~ l3 dS. (3) 
We maximize the log 



(4) 



where the distribution of our measured values of flux and error 
(F/,oy) is given by 



J °°P(F h <j i \S)P(S)dS 
Jf°° J OC P(F i ,a i \S)P(S)dSdF' 



(5) 



A more detailed account of this method, and its variants, can 
be found in Temple et al. (2005). 

Assuming the errors in measuring flux and luminosity are 
distributed as a Gaussian, and that the minimum measured 
flux Fmin is the flux of the faintest source (as quoted above) in 
each observation, we numerically find the value of (3 (thus a) 
for which C in is maximum, individually for both obser- 
vations. 

The ACIS-I data (red histogram in Figure ITOb gives the 
slope of the cumulative luminosity function a = 1 . 14 ± 0.03, 
and the ACIS-S data (blue histogram) yields a = 1.19 ± 0.03, 
which are consistent with each other within errors. From sam- 
ples of 4 and 14 early-type galaxies Gilfanov (2004) and Kim 
& Fabbiano (2004), respectively, found point source luminos- 



Chandra Observations of NGC 4636 



7 



3 

&5 



3 

O 




X— ray Luminosity (erg s ) 



FIG. 10. — Luminosity function of X-ray point sources in NGC 4636 detected in the 0.5 - 2 keV band. Only sourc es w ith > 10 net counts and radial distance 
1.5' < d < 6' from the center of the galaxy are included. The ACIS-I data (above, red histogram) are best fit (see j) |3.4> by a power-law with a = 1.14 ± 0.03 
(dashed line) and the ACIS-S data (below, blue histogram) are best fit by a power-law with a = 1.19 ± 0.03 (dot-dashed line). There are no ULXs in these 
observations. 



ity function slopes in the range a =0.6-1.2; our values for 
NGC 4636 are at the higher end of this range. 
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FIG. 11. — The X-ray point sources detected in NGC 4636, shown by '+'s, 
plotted along with all globular cluster candidates (discrete optical sources 
from Dirsch et al. (2005), brighter than R < 23.5 with C-R color between 
and 2.5), plotted as open circles. The optical extent of the galaxy is shown 
as the D25 ellipse (major and minor axis of 6.0' and 4.7' respectively). Other 
known optical sources (from NED) in the field are also plotted: quasars as 
squares, and background galaxies and clusters as triangles. The supernova 
remnant SN1939A is plotted as a diamond within the optical extent of the 
galaxy. 



3.5. Correlation with the globular cluster population 



It is well-known that a significant fraction of LMXBs in the 
Milky Way and elsewhere is associated with globular clus- 
ters. In the Milky Way, -10% of all bright (> 10 36 erg s _I ) 
LMXBs are found in globular clusters (GCs), even though 
GCs account for < 10~ 3 of the stellar mass of the galaxy (e.g. 
Katz 1975, Clark 1975, Grindlay 1993). This has led to the 
suggestion that LMXBs are very close binary systems formed 
as a result of dissipative two-body or three-body encounters. 
This could be in the form of encounters between neutron stars 
and ordinary stars, which are more likely in the dense cluster 
cores (Fabian et al. 1975; Hut et al. 1992), or the exchange of 
a companion, where a compact object replaces a member of 
a binary in a three-body interaction (Clark 1975, Hills 1976). 
Indeed, it has been suggested that LMXBs are formed primar- 
ily in the cores of globular clusters, and some of the resulting 
binaries are later ejected from their host clusters (White et al. 
2002). If so, one expects the LMXB population in a galaxy to 
be a good tracer of globular clusters. Alternatively, a signif- 
icant fraction of the LMXBs could have formed in the field, 
possibly as part of the last major star formation episode of the 
galaxy (Irwin 2005, Juett 2005), in which case the populations 
of GCs and LMXBs need not be strongly correlated. 

In early-type galaxies, Chandra observations show that 
there is an association of LMXBs with globular clusters. The 
fraction of LMXBs identified with known GCs varies from 
at least 20% in the Virgo elliptical NGC 4697 (Sarazin et al. 
2001) to up to 70% in NGC 1399, the central galaxy of the 
Fornax cluster (Angelini et al. 2001). Furthermore, it is ob- 
served that the LMXBs in early-type galaxies are >3 times 
more likely to be in the redder globular clusters, for galax- 
ies which exhibit bimodality in color of the GCs, e.g. M87, 
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Jordan et al. 2004; NGC 4472, Kundu et al. 2002; more 
recently in ten others; Kundu et al. 2007, Sivakoff et al. 
2007; Cen A (NGC 5128), Minniti et al. 2004, Jordan et 
al. 2007, Woodley et al. 2008; NGC 3379, Brassington et 
al. 2008. Since the redder GCs are relatively more metal-rich 
(e.g. Kundu & Zepf 2007), this can be a diagnostic of the 
characteristics of the compact object's companion star in the 
LMXB, and its history of formation (Maccarone et al. 2004; 
Ivanova 2006). However, this could also result from system- 
atic variation in the initial mass function of the GCs (Grind- 
lay 1987), a possibility that has not been observationally ex- 
plored, but could give rise to the same effect, independent of 
the nature of the companions. 

We matched our list of X-ray point sources with the list of 
globular cluster candidates, from the study of Dirsch et al. 
(2005), which uses a deep mosaic CCD observation using the 
CTIO Blanco telescope. The photometry is available in the 
Washington C and Kron-Cousins R system. Of the optical 
point sources listed in this work which extend over an area of 
0.25 deg 2 around the galaxy, we chose the sources with mag- 
nitude R < 23.5 and color < C-R < 2.5 as globular cluster 
candidates. These are not spectroscopically confirmed to be 
globular clusters belonging to the galaxy. FigureQT]shows the 
distribution of both X-ray and optical point sources, together 
with other background sources that are known in the field. 
In a more recent paper, the group responsible for the orig- 
inal list has published spectral observations of a small sub- 
set of 200 of the original list of GC candidates (Schuberth 
et al. 2006). Of the sources with magnitude R < 23.5 and 
color < C-R < 2.5, > 80% of the candidates were found to 
be globular clusters from measured redshift and the rest fore- 
ground stars. 

The positions of both X-ray and optical sources are known 
to an accuracy of better than an arcsecond. We matched the 
two lists by taking each X-ray point source, and finding its off- 
set from the nearest globular cluster candidate. We checked 
for systematic translation and rotation between the two lists 
by seeking to maximize the matches for small values of rota- 
tion and translation of all sources, but could not improve upon 
the matching done above. Based on the distribution of offsets, 
we chose to limit the search radius to 1.5". If two GC candi- 
dates fell within this radius, we assigned the X-ray source to 
the nearer one in angular distance. Of the 318 sources in our 
list, 77 were matched to globular cluster candidates. These 
sources are marked "GC" in the last column of Tableland 
are listed separately in Table [3] Of these 77 matched point 
sources, 48 lie within 1 .5' < r < 6 .0' of the center of NGC 4636 
and have more than 10 net counts in either of the ACIS-S 
or combined ACIS-I observations. Based on the number of 
X-ray sources and GC candidates in this annulus, we expect 
only three chance coincidences. No other known background 
sources from NED were matched with X-ray sources within 
10' of the center of the galaxy. 

A C-R color histogram of globular cluster candidates is 
shown in the top panel of Figure fT2l (dashed histogram). As 
in many other early-type galaxies, the distribution of the col- 
ors of candidate globular clusters here is bimodel. Dirsch et 
al. (2005) fit a two-Gaussian model to their color distribu- 
tion, and find that the two peaks are at C-R = 1 .28 and 1 .77, 
with the intervening minimum occurring at C-R =1.5. We 
adopt these values, and here refer to those with C-R > 1.5 
as the "redder" GCs, and the rest, the "bluer" ones. Of the 
48 matched GCs between 1.5' - 6.0' of the galactic center, 



13 have colors between 0.85 < C-R < 1.5, and 30 have col- 
ors C-R > 1 .5, so 70% of the matched GCs in Figure [HI with 
colors between 0.8 and 2.5 are in the redder category. 

The difference in color between the two populations is pre- 
dominantly due to a difference in metal abundance, since for 
stellar populations that are more than a few Gyr old, the opti- 
cal colors would be more sensitive to the [Fe/H] index than 
to age (e.g., Worthey 1994, Bruzual & Chariot 1993). In 
NGC 4636, even though bluer GCs are more abundant, a ma- 
jority of the X-ray point sources (LMXBs) are associated with 
the redder GCs (those of near-solar abundance), as is shown 
by the solid histogram in the left panel of Figure [12] which 
represents the color distribution of the GCs matched with X- 
ray point sources. A similar association has been reported for 
other nearby early-type galaxies (Kundu et al. 2002, Jordan 
et al.2004, Kim et al. 2005, Kundu et al. 2007, Sivakoff et 
al. 2007). This is also consistent with the observation that 
most LMXBs associated with GCs in the Galaxy and M31 lie 
in those systems with a near-solar abundance (Grindlay 1993, 
Bellazzini et al. 1995, Bregman 2006). 

Other studies have found that X-ray point sources are pref- 
erentially found in optically more luminous globular clusters 
(Angelini et al. 2001, Sarazin et al. 2003, Kundu et al. 2002, 
Jordan et al.2004, Kim et al. 2005, Xu et al. 2005). GCs 
containing LMXBs also have been noted to be significantly 
denser in the Galaxy and in M31 (e.g. Bellazzini et al. 1995). 
Since the mean size of a GC does not significantly vary with 
luminosity (McLaughlin 2000), more luminous GCs can be 
expected to be denser, and such a correlation can result from 
the higher density of potential companion stars in such sys- 
tems. In the lower panel of Figure [12J we show the distri- 
bution of all GC candidates 1.5' < r < 6' from the center of 
NGC 4636 (dashed histogram), plotted along with those that 
match with X-ray point sources (solid histogram). The dotted 
line shows the ratio of matched to total GC candidates as a 
function of R magnitude. This ratio declines towards fainter 
GCs, indicating that the probability of a GC containing an 
LMXB is proportional to its luminosity, consistent with other 
studies (e.g. Sarazin et al.2003, Jordan et al. 2004, Sivakoff et 
al. 2007, Woodley et al. 2008). 

In Figure [13] we plot the apparent magnitude and color of 
the host globular clusters against the X-ray luminosity of the 
matched X-ray point sources. These plots indicate that the X- 
ray luminosity of the matched point sources does not depend 
on the color or absolute luminosity of the host GCs, which 
agrees with earlier work noted above. Together with the pre- 
vious observation, this should be a useful constraint on the 
formation process of LMXBs in globular clusters. 

Finally, we compare the luminosity function (LF) of the X- 
ray point sources matched with globular cluster candidates to 
the LF of those that are not matched with any optical source. 
Here, we study only the point sources within 1.5' < r < 6.0' 
of the center of NGC 4636, where the X-ray sources are not 
seriously contaminated by background AGN (see Figure 0). 
In Figure [14] we plot the LFs of GC-matched and non-GC 
matched X-ray point sources. The two plots represent X-ray 
luminosities obtained from two different epochs of observa- 
tions, to account for the possible change in LF due to long 
term variability in these sources. The solid histogram shows 
the LF of the matched point sources, and the dashed his- 
togram is the LF of all X-ray point sources that did not match 
a GC candidate. Note that the brightest X-ray sources are not 
matched with any GC candidate. In both cases, the slopes of 
the LFs of the GC sources and non-GC sources in the range 
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FIG. 12. — (Top:) A histogram of the C—R colors of all globular cluster candidates (discrete optical sources from Dirsch et al. (2005) with R < 23.5 and 
< C—R < 2.5), found within 1.5' < r< 6.' of the center of NGC 4636 (dashed histogram), plotted along with those that match with X-ray point sources detected 
in the Chandra observations (solid histogram). It is clear that the X-ray point sources are preferentially associated with the redder globular clusters. (Bottom:) 
Histograms for the same two samples (dashed- all candidates, solid- those that match point sources) as a function of GC magnitude. The red dotted line is the 
number of matched GCs divided by the total number of GC candidates in each luminosity bin. This fraction declines towards fainter GC magnitude. 



1 .8 x 10 37 erg s" 1 < L x < 1 x 10 38 erg s" 1 match within 2a. For 
the ACIS-S observation, from a Kolmogorov-Smirnov (K-S) 
test, the probability that the two samples are drawn from the 
same population is 78%, whereas for the ACIS-I samples, the 
corresponding value is 50%, with the maximum deviation oc- 
curring at the faint end in both cases. We note that the K-S test 
underestimates this probability for discrete distributions (see, 
for example, Sheskin 2003). At the faint end of the luminos- 
ity function, there is emerging evidence, from nearby galax- 
ies like Cen A (NGC 5128) where more complete samples of 
fainter sources can be studied, that there are fewer X-ray faint 
sources found in GCs than X-ray bright ones (e.g. Woodley et 
al. 2008). We see perhaps some evidence for a dearth of GC 
matches at lower luminosities in the ACIS-I LFs (Figure [141 
bottom panel). 



3.6. Individual Sources 
3.6. 1 . Galactic Nucleus 

We detect a soft, luminous (Lx = 2.2 x 10 38 erg s" 1 ) X-ray 
point source at the galaxy's nucleus, as determined from the 
VLA 5 GHz (Griffith et al. 1995) radio position (a(J2000) 
= 190.7077, <5(J2000) = 2.6878). Loewenstein et al. (2001) 
do not detect nuclear activity in NGC 4636 in the 2-10 keV 
band. However, this source is very soft, and ~ 90% of its 
emission falls below the 2-10 keV band. This source is 
shown with a purple X in the color-color diagram (Figure |4]i 
and is significantly softer than the LMXBs. The spectrum is 
well-fit by a power-law model with a photon index of A = 
2.36 (68% confidence range is 1.94 - 2.76). In fellow Virgo 
elliptical galaxies NGC 4472 and NGC 4649, Soldatenkov et 
al. (2003) detect two soft x-ray sources at the galactic nuclei 
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FIG. 13. — The mean X-ray luminosity of the point sources with > 10 net counts that are found to match globular cluster candidates within 1.5' - 6' of 
the galactic center, plotted against the R magnitude (left) and the C — R color (right) of the matched globular cluster candidate. No trend is observed with the 
magnitude of the host globular cluster, while, as seen in Figure [T2l the overabundance of matched X-ray sources in the redder GCs is easily apparent. 



with similar power law indices (2.5 ± 0.4 and < 2.2). 

3.6.2. Additional Black Hole Candidates 

We detected two soft, luminous (L x « 2 x 10 38 erg s" 1 ) 
sources very near the center of the galaxy (d < 5"), whose 
soft spectra, like the source at the nucleus, may be signa- 
tures of black holes. These sources are shown by the pur- 
ple asterisk and plus sign in the color-color diagram (Figure 
EQ. McClintock & Remillard (2004) find that black hole bi- 
naries in the thermal-dominated state are well-described by 
a disk-blackbody model with T = 1 keV. The spectra of the 
two sources are very similar, so we added them before fitting 
with XSPEC (Arnaud 1996). With a disk-blackbody model, 
we find a best-fit temperature of 0.30 keV with a 68% con- 
fidence range of 0.27 - 0.33 keV. The co-added spectra are 
equally well-fit by a power-law model with a photon index of 
A = 2.81 (68% confidence range 2.58 - 3.03). 

Although we are prevented by the insufficient source counts 
from reaching a definitive conclusion, it is possible that these 
sources, and the one at the nucleus, are indeed signatures 
of black holes at the center of NGC 4636. If so, they may 
be massive black holes which fell to the galaxy center un- 
der dynamical friction, in which case they would have masses 
greater than 10 5 M Q (Tremaine & Ostriker 1999). However, 
given the luminosities of these sources, it is more likely that 
they are remnants of merged galaxies, similar to what Ko- 
mossa et al. (2003) observe in the elliptical galaxy NGC 6240. 

3.6.3. SN1939A 

SN 1939A is a Type la supernova 0.57' from the center of 
NGC 4636 (a(J2000)=190.700°, <5(J2000)=2.693° (Zwicky 
1939, Giclas 1939, Tsvetkov & Bartunov 1993). We do not 



detect X-ray emission at this location. By generating a large 
number of Poisson realizations from the observed background 
at this location in our Chandra observations, we calculate a 
3cr X-ray count rate upper limit of 2.5 x 10~ 4 cts s" 1 in the 0.5 
- 2 keV band. This translates to a luminosity upper limit of 
~ 4.5 x 10 37 erg s" 1 in the same energy band. 

3.6.4. Bursting X-ray Source 

A light curve for the bursting X-ray source, detected at 
a(72000) = 190.8825°, (5(72000) = 2.6190°, is shown in Fig- 
ure |9] The matching source (0.8" away) in the 2MASS All- 
Sky Catalog of Point Sources, 2MASS J12433175+0237080 
(a(72000) = 190.882303°, 5(72000) = 2.618898°) has J, H, 
and K band magnitudes of 10.139 ± 0.032, 9.485 ± 0.03, and 
9.293 ± 0.028, respectively (Skrutskie 2006). From the Sloan 
Digital Sky Survey Data Release 6, we get g and r band mag- 
nitudes of 13.92 ± 0.01 and 14.83 ± 0.01, respectively, for 
matching source SDSS J12433 1.75+023708.0 (a(72000) = 
190.8823016°, (5(72000) = 2.61889293°, separation=0.8") 
(Adelman-McCarthy 2008). Based on the source's optical and 
infrared brightness and its large distance from the center of the 
galaxy (1 1.3'), it is most likely not a member of NGC 4636, 
but is probably a relatively nearby flare star in our Galaxy. 

Figure Q3] shows hard and soft band lightcurves for this 
source from the ACIS-S observation. No hard (2-6 keV) 
counts are detected above the background prior to the burst. 

4. SUMMARY AND CONCLUSIONS 

We have analyzed three Chandra ACIS observations, taken 
over three years and totaling ^193 ks, of the Virgo galaxy 
NGC 4636, and have detected 277 X-ray point sources above 
a luminosity of 1.2 x 10 37 erg s" 1 in the 0.5 - 2 keV band, 
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FIG. 14. — The luminosity functions (LF) of all X-ray point sources, found within 1.5'<r<6.0' of the center of NGC 4636 with > 10 net counts that are matched 
with globular cluster candidates (solid histogram). Also plotted are the LFs of all X-ray point sources that did not match a GC candidate (dashed histogram). 
The two plots are for 0.5-2 keV band luminosities calculated for the ACIS-S Sequence 600083 observation (top) and the combined ACIS-I Sequence 600300 and 
600331 observations (bottom). In both cases the LF slopes for the matched and unmatched sources are similar between 1.8 X 10 37 erg s _1 < L x < 1 X 10 38 erg 
s -1 (log=37. 25-38). Note that the brightest X-ray sources are not matched with globular clusters. 



outside the central 1.5' bright galaxy core. Based on the es- 
timated density of AGN in the field, ~ 123 of these are likely 
members of the galaxy, while the rest are likely AGN. In the 
region from 1.5' to 6' from the center of the galaxy, there are 
129 X-ray point sources detected with > 10 net counts, ~25% 
of which are likely AGN. 

We calculate X-ray colors from fluxes (rather than from 
counts, as is traditionally done), and find that this results in 
more clearly grouped populations on the X-ray color-color 



diagram. We identify a large group of LMXBs, and a small 
group of much softer sources, including three within 5" of the 
galactic center which may be black holes. 

We find 77 matches between our X-ray point sources and 
potential globular cluster (GC) candidates found in deep opti- 
cal images of NGC 4636. Choosing the subset of 48 matched 
point sources with > 10 net counts that correlate with GC can- 
didates and lie within 1 .5' < r < 6' of the center of the galaxy 
(out of 129 X-ray point sources in this annulus), we find that 
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FIG. 15. — Soft band (0.3 - 2 keV, top panel) and hard band (2-6 keV, bottom panel) lightcurves from the ACIS-S observation for the bursting X-ray source 
at a(/2000) = 190.8825° , 5(72000) = 2.6190° . Background counts are shown by the dotted line in each panel. 



the overwhelming majority are associated with the redder GC 
candidates, those that are thought to have near-solar metal 
abundance. We find that the fraction of GC candidates with X- 
ray point source matches decreases with decreasing GC lumi- 
nosity. The luminosity functions of the point sources matched 
with GCs and of those that are unmatched have similar slopes. 
We searched for variable sources on timescales ranging from 
hours to years. We find that 54 sources in the common field of 
view (24%) show long-term variability and vary by > 3cr be- 
tween the January 2000 ACIS-S observation and the February 
2003 ACIS-I observations, while about 10% of sources show 
short-term variability within a single observation. 
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9.92 


2 


190.7951 


2.8287 


79 


9.609E+37 


0.06 
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a Net counts and luminosity (erg s ) are in the 0.5 - 2 keV band. 

b Soft color is defined as (M-S)/(S+M+H) where S, M, and H are the fluxes in bands 0.5 - 1, 1 - 2, and 2-8 keV, respectively. 
c Hard color is defined as (H-M)/(S+M+H) where S, M, and H are the fluxes in bands 0.5 -1,1-2, and 2-8 keV, respectively. 

1 means the source was detected by the ACIS-S, 2 means the source was detected by the ACIS-I, and 3 means the source was detected by both instruments. V 
denotes a long-term (i.e. between observations) variable and SV denotes a short-term (i.e. during a single observation, either the ACIS-S (5) or coadded ACIS-I 
(/)) variable. MV indicates that the source flux varies significantly between the sequential ACIS-I observations. GC indicates that the sources is matched with 
a globular cluster. The number following GC gives the source's position in Tablef3] 
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TABLE 3 

X-RAY POINT SOURCES MATCHED WITH GLOBULAR CLUSTERS CANDIDATES IN NGC 

4636 
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0.14 
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0.22 
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1.36 


23.23 
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N 


R.A. (X-ray) 


Dec (X-ray) 


Offset(") 


CTIO ID a 


R.A. (optical) 


Dec (optical) 


C-R color 


R magnitude 


76 


190.6942 


2.7729 


0.70 


6052 


190.69403 


2.77278 


1.73 


21.57 


77 


190.7134 


2.6774 


1.01 


15522 


190.71320 


2.67719 


0.92 


23.01 



a Dirsch et al. 2005 



